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SUMMARY 



This paper presents the resnlts of tests conducted 
to aetermine the effect of the constructional elements 
of a Laval nozzle on the velocity a,ns^, pressure distribu- 
tion and the magnitude of the reaction force of the jet. 
The effect was studied of the shapes of the entrance 
section of the noz^slo and three types of divergent scc^ 
tions: namely, straight cone, conoidal '7ith cylindrical 
end piece and diffuscr obtained computationally by a 
graphical method due to Professor F. 1. Frankly Tho ef- 
fect of tho divergence angle of the nozzle on the jot re- 
action T7as also invest igated^ The results of the investi- 
gation shoTTed that the shape of the generator of the inner 
surface of the entrance part of the nozzle essentially has 
no effect on the character of the flow and on the reaction/ 
The nozzle that T^as obtained by graphical computation assured 
tho possibility of obtaining a flow for which tha velocity 
of all the gas particles is parallel to the axis of symmetry 
•of the nozzle, the reaction being on the average 2 to 3 per- 
cent greater than for the usual conical nozzle under the 
same conditions. For the conical nozzle the maximum reaction 
was obtained for a cone angle of 25^ to 27*^. At tho end of 
this paper a sample computation is given by the graphical 
method, Tho tests were started at the beginning of 1936 and 
this paper was written at the same time. 



IKTRODUCTIOK 



The operation of certain machines and apparatus, for 
example, steam or gas turbines, rocket motors, and wind 

♦Report No. 478 of the Central Aer o-Hydrodynamical Inst''- 
tute , Moscow, 1940. 
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tunnels, is based on the transformation 6f the potential 
energy of the gas contained in a certaia limited volume 
into kinetic energy of the jet flowing out of a container 
with a large velocity through nozzles of one type or 
another • 

For ol:taining jets with su"bsonic velocities cylin- 
drical nozzles are used while conical nozzles (Laval 
noazles) are used for supersonic velocities. 

The theoretical velocity of outflovr is generally 
denoted as tht-.t velocity w' which would be obtained 
if the flow process v/ere a rovnrsiblc adlabatic or 
isentropic process, that is, if the condition were sat- 
isfied 

ds = 0 

or 



Tds = 




P 

v/here 

T absolute temperature of the gas 

6 heat contdnt or enthalpy 

8 entropy 

p pressure 

p density 

Actually the flow process is not reversible and 
isentropic. The transfer of heat to the nozzle walls 
and the friction of tte gas particles at the ^-'alls lead 
to an increase in the entropy so that always ds > 0. 
The true velocity w will, therefore, be less than th<5 
theoretical velocity v/ * • The ratio of these tvo veloc- 
ities is denoted as the velocity coefficient cp = w/w » . 

It should be noted that the transformation in the 
nozzle of the potential energy into kinetic occurs v/ith 
a high degree of efficiency* The velocity coefficient 
cp even for a nozzle v/ith rough inner surface is still 
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near unity and does not drop below 0.92. Tor nozzles 
with a highly polished surface and of a correctly designed 
shape cp reaches a value of the order of 0.98. 

Another important criterion for judging the qvality 
of a nozzle is the ratio of the true discharge 5 in 
unit time to the theoretical discharge G ' . This ratio 
is denoted as the discharge coefficient through the nozzle 

The coefficients cp and p. depend on the initial 
state of thp gas, th^ degree of expansion of the gas in 
the nozzle, and appar^^ntly, on the geometric character- 
istics of the latter. The values of tp and \x for 
each nozzle are doterriiined experimentally. In particular, 
the coefficient cp raay be cOiiiputed as the ratio of the 
momontums obtained exper inontally to the theoretical 
(assunin^^ = 1 ) • C^his assumption is admissible only 
for c oupar is on- purposes . In general, however, \i is 
alv/ays less than unity.) Thus, knowing the discharge 
rate per r^ccond. & , and the force of reaction obtained 
exper iiuentally , and remembering that 



Gw 

P = — = mw 
S 

while the thoorotical reaction of the jot is 



Gw » 

» P » = = mw 



S 

there is obtained 



<p = ^ (1) 



A correctly chosen shape of nozzle should assure 
a flovr with uniformly distributed velocities and static 
pressures both in the transverse and longitudinal sec- 
tions of the flow. The direction of the velocity of all 
the particles of gas should be parallel to the axis of 
symmetry of the nozzle. 

In the case of nozzles for turbines, and rockets the 
satisfaction of these requirements is necessary in order 
to avoid losses in momentum. For nozzles of wind tunnels, 
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however, the condition of parallel velocities of the gas 
particles must "be satisfied since otherwise the character 
of the flow about the bodies tested in the tunnels would 
not he similar to the flow in free flight. 

In this paper the results are presented of tests 
undertaken with the special purpose of Aeternining the 
effect of the constructional elements of the nozzle on 
the magnitude of the rp.actlon of the. jet- and on the char- 
actef of the velocity and pressure dis tr ihut ioji in the 
flow. 



TEST PHOCEDims 

The reaction of the Jet can he measured hy two 
methods. By one method the irapulse of the jet is measured 
\^ith the aid af a disk fliounted normal to the direction of* 
the jet and connected v;ith a dynamometer (fig. l)» Know- 
ing the rate of discharge per second the true velocity 
of the flow can he determined by the momentum equation. 

Thus » if P is the thrust measured by the dynamometer, 
m the mass of air flovring per second assumed constant , and 
w the flow velocity then by definition: 



P = mw 

whence 

P 

w = — 

m 



(2) 
(3) 



The other method is based on the direct measurement 

of the jet reaction expressed by equation (2). In general, 

equation (2) is valid only in the case that p — the pres— 

ure of the gas in the outlet cr o s s— s e c t i on of the nozzle — 

and Pq — the pressure of the external medium— are strictly 

equal to each other. The nozzle satisfying the condition 

T> = p will be denoted as the normal or c onput ed. noz z le . 
^ a 0 

In the case of the inequality J)a, f Pq equation (2) must 
be replaced by the expression 

P = mw+ faCPa-Pt^) (4) 

where fa is the area of the outlet section of the noizzle. 
The last term bn the right side of the above equation is 
the static term. 
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Actually a nozzlo rarely oporatos under the computa- 
tional assumed conditions and therefore the investigation 
of the effect of the static term on the value of the thrust 
is of considerable practical interest. In order to measure 
the reaction of the Jot by the second method an apparatus 
was constructed (general view shown on figs. 2 and S)» 
The tube 1 serves to conduct the comptressed air to cylinder 
2 which, with the aid of a flange, is connected to tube 1. 
Cylinder 2 is surrounded on top by a jacket 3 which at it3 
upper end is supported against the ball bearing 6 and can 
rotate about cylinder 2. To cylinder 3 on the two diamet- 
rically opposite sides are attached two tubes 5 and 5a. To 
tube 5 is attached the weight 9 as a counterweight of chamber 
7 with the tested nozzle 8 attached to tube 5a. A long thin 
manometer tube is attached to connecting tube 11. In the 
case where the working gas consists of the products of com- 
bustion the spark plugs 14 serve to ignite the fuel mixture. 
For the same purpose an opening is provided for mounting a 
spark plug in chamber 7. To make cylinder 3 airtight a 
labyrinth packing is placed between the outer surface of 
the cylinder and the inner surface of the jacket 3. More- 
over, for the same purpose the water between cylinders 2 
and 3 was put under a certain pressure, greater than the 
pressure in cylinder 2. The water is supplied through tube 
13. In working with the products of combustion this water 
serves also to cool cylinder 2. Practical operation of the 
apparatus sho-^ed that the hermetic scaling of cylinder 2 
was assured by the labyrinth packing alone and there was no 
neodof a water back pressure. The fiber plate 10 serves as 
a safety valve in case of a sudden rise in pressure in cyl- 
inder 2. Chamber 7 is supported at the bottom on the dyna- 
mometer with the aid of which the reaction force is measured. 
The apparatus also permits the carrying out of a qualitative 
comparison of two nozzles. Per this purpose, in place of the 
plug 12 in chamber 7, a second nozzle is screwed in the re- 
action of which is directed opposite to that of nozzlo 8. 

To measure the rate of discharge of the air through 
the nozzle a so— called diaphragm is placed in pipe 1 (not 
shown on the sketch). The pressure difference h at the 
inlet and outlet of the diaphram is measured by a differ- 
ential manometer of the Yasinsky type. The first term In 
expression (4) for the reaction force of the jet, that is, 
the momentum^, is measured n^ar the nozzle by a dynamometer 
with disk placed normal to the nozzle axis, the distance 
from the outlet section of the nozzle to the center of the 
disk should not be greater than two diameters of the jet 
in order to avoid the unfavorable distorting action of the 
disturbed outside air sucked in by the jet. 
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The dynamic head was measured with the aid of a 
sytstem of Pitot tubes mounted on the disk for measuring 
the momentum (fig. 4)^ the distance between the individual 
tubes bein,^ equal to 5 millimeters. By displacing: the 
tubes along the axis of the jet the true dynamic pressures 
along the diameter and along the axis of the jet may be 
determ ined. 

Knowing the value of the dynamic pressure for a given 
point it is not difficult by the known formula of Raleigh 
(reference 1) to obtain the true velocity. From numerous 
tests it was found that the values of the velocities ob- 
tained by this formula are very nearly the true flow 
velocit ie s, 

Figure 5 gives the curve, computed by the Raleigh 
formula, of the dynamic pressure in a supersonic jet flow 
as a function of the Bairstow (or Mach) number (w/c). 
This curve may be used to determine the actual velocity 
of the supersonic flow (or the value of w/c) from the 
dynamic pressure experimentally obtained. 

The temperature in the chamber ahead of the diverging' 
nozzle was measured with a thermocouple. The temperature 
of the air in the flow was determined on the assumption of 
adiabatic expansion of the gas in the nozzle by the formula; 

k-X 

T = Ti ^i) (7) 

where 

T^ and p^ temperature and pressure of the gas in the chamber 

T and P corresponding parameters after expansion 

k ratio of specific heats 

The static pressure in the flow was measured by a tube 
placed in the gas flow with its axis parallel to the velocity 
direction. Figure 6 shows the static tube mounted in the flow. 

Besides determining the above thermodynamic parameters 
the flow spectrum was also phot o^rraphed by the Tepler method 
during the tests^ 
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&AS JET OF SUPERSONIC VELOCITY 



In investigating the motion of gases with small 
velocity the phenomenon of comprGssihility is generally 
not taken into account. The error in the pressure com- 
putations of the flow with the compressibility neglected 
constitutes ahout 3 percent at speeds of 100 meters per 
second. At speeds of the order of 200 meters per second, 
horrever, this error is more appreciable and constitutes 
about 9 percent while for w 340 meters per second, 
that is, for flow velocities equal to the velocity of 
sound, the error reaches 27.5 percent. 

Figure 7 shows the lavr of velocity distribution in 
the transverse cr os s— sc ctt i on and figure 8 the velocity 
and static pressure distribution alona; the flow axis 
for a conical Laval nozzle. As may be seen from the 
curves the velocities and static pressures in both in 
the transverse and axial sections are not constant^ The 
sharp, almost periodical fluctuations in the values of 
the static pressures arc explained by the fact that in 
the sonic and supersonic flow a very complicated system 
of stationary -raves arises. Those waves -rere first ob- 
served by Ernst Mach in whose honor they are called Mach 
waves. 

To explain the physical nature of the Mach waves 
a two— dimensi onal flow with velocity greater than that 
of sound is considered (fig. 9). It is asa.umed that the 
pressure in the outlet cross section of the nozzle p^ 
•is greater than the pressure Pq of the external medium. 
At the outlet section at point A a disturbance then 
arises as a result of the discontinuous pressure drop from 
Pa Pq» "^-^^ rarefaction wave having a velocity always 

equal to the local velocity of sound is propagated only as 
far as the line AB, The latter is a Mach wave and is essen- 
tially the envelope of all the sound -raves arising at point 
A. The Mach ^ave AB is thus the boundary of the rare- 
faction disturbance starting at point A. The expansion of 
the gas in the region ABC frOm p^^ to p^ produces a 
change in the radial velocities of the florr particles as 
a result of -^hich the direction of the flo^ from line AB 
deviates toward the outside until the pressure is equal to 
that of the outside medium (line AC) after which the flow 
again becomes straight. In the case -rherc the pressure in 
the outlet section of the nozzle is less than the pressure 
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of the outside medium, that is, p^^ < p^, the deviation 

of the flow will he directed inward,. Therefore in the 
first case, that is, Pa Po diameter of the flow 

will bo groater than the diameter of the outlet section 
of the no^i'^le while in the second case, p^^ < p^ the 

diamete:: the flow will be less than that of the out- 
let secti:i-i In this case the line AC is absent from 
the flow 3p^^r:trum. T>.us in a flow with axial symmetry 
the expansion due to the inequality of p^^ and p^ 
will occur ill the region included between the two conical 
Mach surfaces with the corresponding generators AB and 
AC, In the region lying ahead of the Mach line AB the 
disturbing effect of point A is entirely absent and the 
pressure in this part of the flow remains equal to the 
pressure in the nozzle, cross-— secti on,, 

The anglo of inj?.inaticn of the Mach wave a is 
determined by ':he ratio of the local velocity of sound 
to the velocity of tha gas before expansion and the angle 
a' by the same ratio after expansion, that is, 

sin a - ' (8) 

w 

The Mach wave diagrams of the supersonic plane- 
parallel flow for Pg^ > Pq and p^ < Pq ^i^^ shown in 
figure 10 and 11. The spectrum of the flow for p^^ > p^ 

is considered. The rarefaction wave AC on reaching the 
boundary layer of the jet is reflected toward the inside 
of the jet in the form of compression waves (the compres- 
sion waves are drawn full line and the rarefaction waves 
dotted)^ In regions 2, 3, 5, and 6 the static pressure is 
equal to the pressure of the surrounding medium, while in 
region 4, separated from the outlet section of the nozzle 
by a distance of two rarefaction waves, the prevailing 
pressure will be lower. In region 7, immediately behind 
the compression waves, the pressure will bo the same as 
in region 1, and there is, therefore, a definite disturb- 
ance period present. Analogous considerations apply to 
the relations of figure 11. The waves AC and A^C^j^ in 
this case will be compression waves and there will he 
corresponding regions of increased and lowered pressures. 

The above description of plane— parallel flow cannot 
without reservations be entirely identified '^ith the flow 



NACA TM No. 1066 



9 



with axial syninjetry because the actual floT7 is three 
dimensional and a Mach wave will bo formed at each point 
on the circumference of the outlet section. 

The Mach waves may be formed not only in the flow 
behind the outlet section of the nozzle but also w.thin 
the nozzle. The latter waves will occur only in case 
the nozzle has a large divergence angle, as a result of 
which the flow separates from the walls even within the 
nozzles, or in case of a very rough inner nozzle surface. 

On the basis of what has already been said the con- 
clusion may be drawn that the nonunif ormi ty of the static 
pressure distribution in the supersonic flow is explained 
by the presence of Mach waves in the flow. This conclusion 
is entirely confirmed also by the results of the tests con- 
ducted. Besides the effect of the stationary sound waves 
the friction of the gas particles at the nozzle walls and 
at the external medium also affects the transverse and 
axial velocity distribution. 

There is still another factor that must be mentioned 
which affects the velocity distribution of the flow^ A« 
test results have shown, the diameter of the flow, with 
increasing distance from the outlet section of the nozzle 
increases as a result of the sucking in of the outer air 
by the jet. This suction takes place with simultaneous 
decrease in velocity along the jet. To investigate this 
phenomenon the following tests wore conductod. A screen 
of thin constantan wires stretched on two porcelain rods 
was mounted near the jet escaping from the Laval nozzle. 
The wires were heated to a red heat by an electric current : 
passed through them. The Teplcr photograph (fig. 12) 
shows clearly the direction of motion of those air particles 
which before being sucked into the stream passed/by the 
heated wires. As may be seen from the photograph the suction 
begins at the outlet section of the nozzle. TSviiently, with 
increasing distance from the nozzle tho quantity of air sucked 
Into the air stream increases and therefore the flow velocity 
correspondingly decreases, the decrease being greatest at the 
boundary layer. 

It should also be mentioned that the flow obtained for 
nozzles with p^^ 4 Pq ^iH always have nonuniform velocity 
and pressure distributions as a result of vortex formation 
and the dissipation of the jet at the outlet section of the 
nozzle. 
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EFFECT OF STRUCTURAL BLBkENTS OF TH^S NOZZLE ON 
TgE THRUST AND CHARACATi^R OF THE FLOW 



The nozzle cavity (fie:, 13) is divided into two 
parts. The first part converges toward the narrowest 
critical section while the second part is a diverging 
cone having an 8^ to 12^ cone angle. The effect of each 
of these parts on the character of the flow and on the 
magnitude of the reaction force of the jet with the ob- 
ject of finding the most rational shape of nozzle is to 
he considered. 



Converging Part of Nozzle 

In this case if the flow into the nozzle is from a 
vessel whose volume is many times greater than the volume 
discharged per second the velocity of the gas at the con- 
vergent section is generally neglected so that the shape 
of the first part is likewise not taken into account. The 
only requirement for this part of the nozzle is that it 
should he carefully rounded and that it should assure a 
smooth flow of the gases. The smoothing curve is chosen 
on an intuitive basis. Considering, however, the motion 
of the gas in this part of the nozzle as analogous to the 
motion of a fluid in a curved channel it may he concluded 
that as a result of the curvature of the path centrifugal 
pressures occur in the flow. The direction of these pres- 
sures will always he a-^ay from the center of curvature. 
Hence in the first part of the nozzle the pressure near 
the axis will ho greater than at the walls. Near the crit- 
ical section the opposite will be t ruo since at this place 
the curvature of the walls is reversed with the center of 
curvature outside the nozzle cavity. It is thus clear that 
the pressure and velocity distributions in this part of the 
nozzle will depend on the radii of curvature of the inner 
surface of the nozzle. From invc s t i -rat i ons of the nozzles 
of wind tunnels of supersonic velocities in various aero- 
dynamic laboratories the dependence of the character of the 
flow on the shape of tho generator of the inner surface at 
the entrance part of the nozzle was cstablishodt Tho nozzles 
of wind tunnels of subsonic velocities have no second (di- 
ve rgent ) part . 

Figure 14 shows the distribution of the dynamic pres^ 
sures q in the flow cross— sect ion for two noszles with 
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different shapes of entrance parts and fieuve 15 sho'vs the 
distribution of the dynamic pressures q and static pres- 
sures X along the axis for the same nozzlofu (reference 2)* 
The dotted curves ^ere obtained for nozzles Tor which the 
generator of the entrance part is a lemniscate and at the 
outlet is an arc of a circle (CAHI T— 5)* The continuous 
curves were- obtained for nozzles whose inner contours were 
constructed by the formula 




where 



Pq radius of the outlet section 

radius of the inlet section 

r variable radius at the distance z from the origin 

of coordinates 

a/ s/^ length of the projection of the curve on the z— axis 

Prom the tests conducted at CAHI it was found that it was 
best to take 

^ = 4.14 ro 

This formula was obtained by litoszynski (reference 3) 
and represents the case of a particular solution of the 
Laplaco equation for the stream function (reference 3). As 
may be seen from figures 14 and 15 the velocity and static 
pressure distributions in the flow are more satisfactory 
for the nozzle computed by litoszynski. 

With the object of determining the effect on the re- 
action force of tho shape of the first part of a nozzle 
having also a diverging part (Laval nozzle) the following 
tests were conducted. Three nozzles were prepared: 1, 2, 
and 3 (fig* 16), for which the diverging parts were entirely 
the same, that is, the diameter of the outlet section da, 
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the diameter of the critical section and the cone 

angle a '^ere equal but the shapes of the entrance parts 
V7ere different^ The shape of the curve for the entrance 
part of nozzle 1 '^as constructed according to the Witoszyn— . 
ski equation. Nozzle. 2 differed from nozzl© 1 in that the 
former had two rings at the entrance part. Finally, nozzle . 
3 had a nonst reamline shape for the entrance part. All 
three nozzles '^ere computed for tt/c = 2^.64. The thrusts 
of the nozzles were measured by the reaction method on the 
apparatus previously described. In table 1 the values of 
the measured thrust and the velocity coefficients cp com-* 
puted by formula (1) are given. 



TABLE 1 



Kozz le 


P kg 




Percent 


1 




0, 965 


100 


3 


29. 9 


.953 


98.7 


? 


29.8 


,950 


98.2 



The magnitude of the theoretical impulse P' was 
determined on the basis of the following: considdrat lons» 
The thermodynamic parameters of the air in the chamber 
during the test wore taken as the following: pj = 21.5 
atmospheres, Ti = 273^. Hememberin:? that the flow is 
open the static pressure in the flow may bo taken approxi- 
mately «qual to the pressure of the external medium, that 
is, p = 1 atmosphere. The temperature of the air in the 
flow was computed by formula (7) 



k~i 




The theoretical absolute velocity of the flow was 
found from the equation 




(10) 
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and the theoretical discharge per second through the 
i^ozzle obtained froiij the expression 

0 = (11) 



T^here v is the specific volume of gas in the flow. Sub- 
stituting the numerical values there is obtaiaed 

T = 115^, w> = 665 meters per second, Gs= 0.545 kilogram 

per second. Hence the theoretical thrust of the notzle is 
equal to 



0 545 

pi = m77» = ~ 555 = 31»4 kilograms 

9*81 



The maximum thrust obtained from the tests on nozzle 
1 T7as equal to 30,3 kilograms. As may be seen from table 
1 the values of the velocity coefficients obtained for 
all throo nozzles do not exceed the limits of th« values 
usually assumed in the computations. For nozzle 3, having 
an entrance part of unfavorable shape the velocity coeffi- 
cient cp is 1.8 percent less than the value 100 assumed 
for nozzle 1; for nozzle 2 cp rras only 1.3 percent less. 

It must be said that the effect of an improvement in 
shape of the first part of the nozzle on the coefficient cp 
is negligible and lies Txithin the limits of experimental 
error. Remembering, hovrever, that the values of cp were 
obtained as the average values of a larce number of tests 
(not le^s than 15 tests for each nozzle) these values may 
be considered as entirely satisfactory. 

The flo^ spectra obtained by the Tepler method are 
similar for all three nozzles, Evidently a chancre in 
shape of the first part of a Laval nozzle has no effect 
on the velocity and static pressure distributions in the 

f lOTT, 



Divergent Part of the Nozzle 

Depending on the ratio of dimensions of the outlet 
and critical cross-sections in the conical part the expan- 
sion of the gase« is accompanied by a drop in pressure from 
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that at the critical section = 0.538 Pi (for p.ir) 

to the outside pressure^ In the case of nozzles for 
which the ratio of the cone diameters ^as assumed e^ireator 
than tho design value the pressure inside the cone may 
drop bel077« the pressure of the external medium (back pros- 
sure As a result of this expansion the pressure rise at 
the outlet occurs abruptly being accompanied by a loss in 
velocity. 

The character of the velocity and static pressure dis- 
tribution in the flo-r from a conical nozzle depends chiefly 
on the structural elements of the second part, that is, on 
the ratio d^^/d^ or, for a given length of nozzle on the 

divergence angle a and also on the contour of the inner 
surface of the nozzle. 

Possession of a general solution of the problem of the 
optimum shape of the diver^^icnt part of the nozzle is not 
yet obtained. Some investigators, it is true, have succeeded 
in solving this problem, experimentally, for certain condi- 
tions, Stanton found empirically the shape of the nozzle 
for the supersonic velocity 'Jind tunnel at the National 
Physical Laboratory, This nozzle gives a flo^r free from 
Mach rraves in the 'working part of the tunnel* The Stanton 
nozzle differed from the usual Laval nozzle in that tho 
divergence angle^of the conical part did not exceed 5^ as 
compared with 10 to 12^ usually assumed for the Laval 
nozzle. In addition tho conical part of tho nozzle smoothly 
passes over into a cylindrical piece at tho outlet. 

Figure 17 shoTrs the curves of variation of the static 
pressure p and velocity vr for this nozzle. As may be 
seen from the curves the cylindrical part has an interval 
whore the curves p and w are almost parallel to the 
nozzle axis, that is, constant. This part was the one 
chosen for the test model. Tho dotted curves on tho fig- 
ure give- the nozzle shape computed under conditions of 
adiabatic expansion of the air in the nozzle without taking 
account of the decrease in velocity due to friction at the 
wall. This confirms the adiabatic law of expansion of the 
air in the nozzle. 

With the object of studying the effect on tho charactor- 
of the flow of the shape of tho curve generating the second 
part of the nozzle a number of nozzles with conoldal inner 
surface and cylindrical piece at tho outlet wore investigated. 
The nozzles were tested for w/ c = 1.41, 1.82, 2.24, and 2.74. 
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The rogults of all the tests ^ere similar and thcrcforo 
the test results for only one nozzle for t/c = 1,82 
arc given. The static pressnres measured along tho floT7 
axis at various heights (H mm) from the outlet section are 
given for this nozzle in table 2, The temperature of the 
air in the flo';7 behind the nozzle was computed by formula 
(7) on the basis of the measured air temperature in the 
chamber , 



TABLE 2 



Test 


H 
(mm ) 


0 

ti c 


0 

t c 


P 

(atm ) 


Test 


H 
(mm ) 


0 

t^ D 


0 

t C 


(atm ) 


1 


5 


190 


8 


0.742 


11 


50 


0 


-110 


1,20 


2 


15 


195 


9 


.749 


12 


60 


0 


-110 


1.40 


3 


25 


185 


3 


1.10 


13 


60 


0 


-110 


1.40 


4 


?o 


200 


13 


.667 


14 


60 


0 


-110 


1.40 


5 


50 


185 


3 


.680 


15 


70 


0 


-110 


1.25 


6 


40 


195 


9 


1.10 


16 


70 


0 


-110 


1,25 


7 


40 


0 


110 


1.20 


17 


80 


0 


-110 


.910 


8 


40 


0 


-110 


1.20 


18 


80 


0 


-110 


. 922 


9 


50 


0 


-110 


1.20 


19 


90 


0 


-110 


.624 


10 


50 


0 


-no 


1.20 


20 


90 


0 


-110 


,624 



From the data of table 3 the curve of static pressure dis-^ 
tribution along the flo^ axis was dra'^n (fig, 18)^ As 
may be seen from the curve the static pressure distribution 
Is extremely variable. The pressure gradient at the center 
of the jet reaches a value of the order of ±0,4 atmosphere^ 
In order to eliminate the effect of possible accidental 
errors the measurements for a ferr point s were repeat ed on 
various days. The. scatter of the test points obtained in 
certain cases may apparently be explained by the effect of 
the variable atmospheric conditions (temperature and baro- 
metric pressure) not taken into account. 

Comparing this curve with that of the static pressure 
curve at the center of the flow obtained by Stanton for tho 
usual conical nozzle (fig. 8) the agreement of the law of 
variation of static pressures along the axis of the flow 
can be readily established. The curve given here is 
reversed (mirror reflection) with respect to the Stanton 
curve, a fact which is evidently explained by the different 
procedure in evaluating the test results. Figure 18 also 
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shows the curve of velocity distribution (w/c) measured 
alon^ the flow axis. As may "be seen from the curve the 
velocity at the center of the jet over a distance from 
50 to 90 millimeters from the outlet section is almost 
constant. For the conical nozzle (fig.. 8) the curve 
w/c cannot be assumed constant over any part of the 
flow. 

Tests were also conducted for measuring the static 
pressure Pj, at various distances r from the flow axis. 

The results of these tests are given in table 3« 



TABLS 3 



H 


r 




P 


(mm ) 


(mm ) 


(atm ) 


(atm ) 








(center ) 


40 


10 


1.00 










\ 1.20 


40 


6 


1.00 


J 


60 


7 


1.26 


■\ 








1 1.40 


60 


10 


1.20 




70 


10 


1.20 


1.25 



The pressures at the center of the jet are greater 
than the pressures at the periphery. 

The cross— sect ional velocity distribution for the 
same nozzle, as may be seen from the experimental test 
curves (fig. 19 ) , is more satisfactory than for the conical 
Laval nozzle (fig, On this basis it may be concluded 

that the conoidal shape with cylindrical end piece gives 
a considerably bet ter (cross— sect ional and axial ). velocity 
distribution than the Laval nozzle. The nonuni f orrai t y of 
the static pressure is, however, also maintained for the 
conoidal nozzle. 

In 1934 Frankl (reference 4) published a method of 
graphical construction of the second part of the nozzle. 
This method is based on the theory of charact eri-st ic s and 
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is an application of the method of Prandtl and Busemann, 
developed "by them for plane— paralle 1 flow, to the case of 
flow with axial symmetry. 

Figure 20 shows a nozzle constructed graphically "by 
the Frankl method. The dotted lines indicate the contours 
of the usual conical nozzle with the same ratio d^^/d^ 

as for the Frankl nozzle. 

Figure 21 shows the curves of the static pressure 
gradient along the flow axis for the Frankl nozzle for 
w/c = 1.47. The pressure gradient for the design condition 
(p^ = 3.49 atm) does not exceed i0.02 atm, . For the case 
where the nozzle operates at other than design conditions 
(pj ^ 4 atm) the axial static pressure distribution is less 
satisfactory as compared with the design conditions. The 
same figure gives the curve of velocity variation along the 
flow axis. As may he seen from the curve the velocity re- 
mains constant at a distance of 30 to 90 millimeters from 
the outlet section. 

For the same nozzle tests were conducted on the trans- 
verse velocity distribution. The results are given on fig- 
ure 22. For the design conditions (full curves) the veloc- 
ities over the cross-section are uniform with a slight drop 
at the boundary of the flow. Near the exit section of the 
nozzle the nonunif ormi t y of the velocities is greater. The 
dotted curves show the velocity distribution for pressures 
of 4, 5, and i atmospheres. As may be seen from the curves^ 
the Frankl nozzle for other than design conditions give^ a 
flow with nonuniform cr oss— se ct i onal velocity distribution. 

! 

For the same nozzle the reaction of the jet was measured 
It was found that the thrust for the nozzle computed accord- 
ing to Frankl is approximately 2.5 to 3 percent greater than 

the thrust of the usual conical nozzle with same d /d 

a c 

ratio. In the case of application of this nozzle to rocket 
apparatus, however, it is necessary to take account of the 
fact that the length of the Frankl nozzle is almost tw^ce 
its equivalent conical nozzle. Therefore the heat losses 
over the nozzle surface during prolonged operation of the 
rocket motor will be greater than for the conical nozzle* 

Simultaneously with the tests for w/c = 1.47 tests 
were also conducted for the Frankl nozzle f#r w/c = 1%47, 
1.62, 1.76, 2.06, and 2.64. The curves of dynamic and static 
pressure distributions obtained for these nozzles are simi ar 
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to the curves for the nosszle for w/c = 1,47, The pressure 
gradient along the flow axis for tt/c = 2,64 is given in 
table 4, 



TAB LIS 4 • 



H (mm ) . 


0 


20 


40 


60 


80 


100 


; P (atm) 


0,95 


0,90 

- 


1,00 


0, 99 


0,90 


1.02 



The cross— sect i onal velocity di s t rilDut i on for the 
same nozzle is shown on figure 2?, Unfortunately^ the 
investigation of the nozzle flow for w/c = 2.64 was 
carried out at the sections near the outlet of the nozzle 
(20, 30, and 40 mm) and, as previously stated, in this 
region of .flow, for the nozzles computed by the Prankl 
method, a certain nonun i f ormi ty is observed in the veloc- 
ity distribution. As the distance from the outlet edge 
increases, however, the character of the velocity distri- 
bution improves* 

The lack of uniformity of the velocities at the out- 
let section of the nozzle may be explained by the fact that 
the back pressure at the outlet is not exactly equal to the 
computation value due to the variation of the barometric 
pressure. For this reason near the outlet section disturb— 
ances ajppear which affect the character of the flow. In 
order to assure the rated operating conditions of the 
Frankl nozzle it is necessary each time to choose the mag- 
nitude pj V as a function of the atmospheric pressure main- 
taining the design ratio p^/p^^* 

As may be seen from the curves of figure 23, on ap- 
proaching the boundaries, the velocities are not lowered 
for the part of the flow investigated and the test points 
lie on on© line. This is evidently explained by the fact 
that at large velocities in the flow region near the nozzle 
outlet section no mixing of the gas of the boundary layer 
with the outside air occurs. The boundary layer particles 
having a large forward velocity form as it were a firm wall 
relative to the outside medium and the flow occurs with 
small friction in the tube formed by the surrounding air. 

In passing from the computation pressure in the chamber 
Pi = 21.5 atmospheres to 22 atmospheres the uniformity of 
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the velocities is maintained ^ith a slight increase in 
the velocity in absolute value* 

Figures 24 and 25 show the flo-TT spectra obtained for 
w/c =s 2,64 with the conical Laval nozzle and with the 
nozzle computed "by Prankl, Both nozzles have the same 
ratio d^/d^ and the same pressure in the chamber p = 

a, C 

21,5 atmospherest In order to maintain the same external 
conditions (atm pressure and air temp.) both photographs 
were taken on the same day. As may be seen from the photo- 
graphs, instead of the complicated system of stationary 
waves formed in the flow of the conical nozzle only weak 
Mach waves are observed in the Trankl nozzle. Moreover 
the flow in a latter remains cylindrical over a considerable 
distance from the nozzle. In the ideal case, that is, for 
stric^t maintenance of the ratio Pq/pi careful polish- 

ing of the walls the Frankl nozzle should not give any Mach 
waves in the flow at all. 

It has already been shown that to satisfy accurately 
these conditions is almost impossible due to the fluctua- 
tions of the atmospherical conditions. Hence the Mach 
waves in a supersonic flow are unavoidable; but their 
intensity, as may be seen on figure 25, is very small in 
the Frapkl nozzle. 

In conclusion it may be said that of the three possible 
shapes considered of the second part of the nozzle; namely, 
(1) the conical, (2) conoidal with cylindrical piece, and 
(3) nozzle cbmputed by Frankl, the latter best satisfies 
the requirements imposed on the nozzle. An example of a 
nozzle computation by the Frankl method is given. 



Effect of the Divergence Angle 

From what has previously been said it is clear that 
the cone angle of the second part of the nozzle may be 
either constant along the nozzle axis or variable. The 
most favorable nozzle from the point of view of the heat 
losses should have the maximum possible divergence since 
the length, and surface area will then be the least. 

With the object of investigating the effect of the 
cone angle on the reaction a number of nozzles with con- 
stant ratio da/dc = 1.73 were tested, the cone angle being 
varied from 10^ to 30^. The pressure and temperature of th 
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air in the cham'ber for all tests 'wrere constant and equal; 
= 21.5 atmospheres, = 273^ absolute. 

Figure 26 shows the dependence of the pozzle thrust 
on the cone angle for rated and nonrated operating con- 
ditions. For the assumed initial thermodynamic parameters 
of the air and for the design conditions the most favorable 
value of the cone angle evidently lies near 25 , since with 
further increase in the angle the thrust curve remains al- 
most parallel to the axis of abscissas. It may he stated 
that the magnitude of the reaction depends very little on 
the divergence of the nozzle. This conclusion agrees with 
the results of the tests conducted under different condi- 
tions, in particular for nons tat i onary flows. To arrive 
at more definite conclusions it would be necessary to con- 
duct tests on nozzles with cone angle greater than 30^ in 
order to determine the point of inflection of the curve of 
reaction force. 

In addition to the tests described above, measurements 
were made on the reaction of the jet for a nozzle with cone 
angle a = 12^ for various values of d^/d^, obtained by 

simply shortening the same nozzle. The results are given 
on figure 37. The theoretical values of the thrust P^j^ 

are computed by formula (4). The maximum test value of 
the thrust was obtained for the design ratio dg^/d^ - 1«73 

and was equal to 18.8 kilograms. For ^^f^c ~ ^» that is, 

in the absence of any divergence the thrust was equal to 
17 kilograms, that is, only 11.0 percent less than for the 
design nozzle. Therefore in working with a cold gas, when 
no heating occurs, the part played by the divergence of the 
nozzle in producing the thrust is relatively small. 

These tests also made it possible to determine the 
character of the changes in the two terms entering the 
expression for the nozzle thrust (equation 4). The results 
are given in figure 28 from which it is seen that with in- 
creasing value of Pq^ the value of the component nw 
rapidly drops. On table 5 the relative values of both 
components for certain values of pg^ are given. 
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TABLE 5 



Pg^ (atm) 


1 


1.6 


2. 8 


5.4 


ll.O 


P (U:g) 


22.7 


22.5 


22.2 


21.0 


20.1 


mw : P 


1 


0. 95 


0,89 


0.80 


0.605 


fa(Pa-Po) = P 


0 


0,05 


0. 11 


0.20 


0,395 



As may "be seen from the tatle the value of the static 
term increases with increasing and for p^^ = 11«2 

atmospheres constitutes 39,5 percent of the total thrust 
of the nozzle. For the usual conditions, however, when 
the lowering of the hack pressure, due to inaccurate com- 
putation of the nozzle or to a change in the atmospheric 
conditions, does not exceed 0,5 to 0.6 atmospheres the 
magnitude of the static. term does not exceed 5 percent of 
the total thrust of the nozzle. The total nozzle thrust 
slowly decreases with increasing Pg^ due to the losses 

in the discontinuous pressure change at the nozzle outlet. 



CONCLUSIONS 



The results of the conducted tests provide a hasis 
for certain conclusions which should he taken into account 
in selecting the shape of the nozzle: 

(1) The shape of the entrance part of the nozzle does 
not show any large effect on the character of the flow and 
therefore tho choice of any particular shape is conditioned 
chiefly hy the structural considerations. 

(2) The nozzle with exit part conoidal and ending with 
a cylindrical piece gives a hetter velocity distrihution in 
the transverse and longitudinal sections of the flow than 
the Laval nozzle hut with regard to the static pressure dis- 
tribution it does not differ from the usual conical nozzle. 

(3) The nozzle computed hy the Frankl method gives a 
parallel flow with uniformly distributed velocities and 
static pressures both in the transverse and longitudinal 
sections. The reaction of the jet obtained with this nozzle 
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is on the average 2 to 3 percent greater than in the usual 
nozzle under the same conditions. 

(4) The maximum reaction is olstained for a cone angle 
of 20^ to 25° and for the normal value of d^^/d^, that is, 

under the condition that the pressure at the outlet section 
of the nozzle is equal to the pressure in the surrounding 
medium. In the ahsence of the second part of the nozzle 
the reaction is found to be 11.0 percent less than for the 
computed nozzle. 

EXAMPLE OF NOZZLE COMPUTATION BY THE FRANZL METHOD 



The theoretical iDasis of the graphical method 6f con- 
structing a nozzle of axial symmetry together with a short 
discussion of its practical applications is given in the 
Prankl article (reference 4). In order to make this method 
available not only to engineers hut also to less qualified 
workers a detailed numerical computation of a nozzle for a 
velocity w = 900 meters per second is given. 

The nozzle (fig# 29) is divided into four parts: 

(1) The part of the nozzle in which the air is con- 
ducted to the critical section. 

(2) The critical section (throat) 

(3) The part between the critical section and the 
streamline part of the nozzle, the generator of the sur- 
face of this part being a straight line (cone) 

(4) Streamline part of the nozzle whose generating 
curve is obtained by the graphical method 

The shape of the curve of the first part of the nozzle 
is not computed but is chosen arbitrarily on the basis of 
considerations of smooth conduction of the air to the critical 
section. To compute the critical section it is necessary to 
have data on the initial condition of the air in the chamber 
(before expansion in the nozzle) and know the rate of dis- 
charge per second through the nozzle. These data are gener- 
ally obtained by a preliminary computation. The initial 
parameters for computing the nozzle are denoted by the follow- 
ing symbols: 
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rato of discharge of the products of combustion, 
kg/ soc 

initial area of the entrance part, 

ressure (pressure in the chamber), kg/cm 





initial 


Pi 


init ial 




initial 




initial 




initial 



With these values given the area of the critical 
section f^ and all the thermodynamic parameters associated 
with it ar'^' obtained; 

Tc • and 



Since in rocket motors and in wind tunnels of large 
velocity it is desirable to have velocities much greater 
than the velocity of sound 

_k_ 

, ^ « / 2 \k-.i 
— > 6 whe re 6 = ( ) 

Pi ^ K^iJ 

With this in mind use is made of the corresponding formulas 
of the flow of gases, given in any book on thermodynamics. 
The process of expansion of the gas in the nozzle is assTzmed 
to be adiabatic. Froa the formula 




(1) 



the area of the critical section of the nozzle may be 
determined. The parameters of the gas in the critical 
section are determined from the following expressions: 



w = y 2g piVi = /2g RT. (2) 

/ k+1 / k+1 
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/ 2 \- 



BIT 




(3) 





(4) 




(5) 



Beyond the critical section the third part of the 
nozzle begins in the form of a diverging cone. The cone 
angle a is generally chosen from 8^ to 13°. The fourth 
part of the nozzle is computed graphically. In order to 
make use of the method of characteristics in constructing 
the boundaries of the flow, it is necessary to know the 
value of the velocity in the direction of the lines of flow 
at the start of the fourth part and the area in this section 
and the direction of the wall of the third part must be 
known. Knowing the cone angle the position of the wall of 
the third part is determined and therefore also the direction 
of the streamlines at the beginning of the fourth part of 
the nozzle. 

Taking any section x at a certain distance from the 
critical section and being given the ratio of the velocity 
w^ in the section x to the velocity in the critical 
section w^, that is, ^x^^;;» area of the nozzle at 

section x is found, that is, at the beginning of the 
fourth part from the following relations: 



f = f 
•^x ^c 




(6) 



or 




(7) 



The entire graphical computation of the nozzle is con- 
ducted with the aid of two graphs: (1) the vector diagram of 
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the velocities in the system of coordinates Vp, v^, 
the radial velocities being laid off on. the axis 
and the axial velocities on the vz axis, and (2) the 
graph of the axial section of the nozzle (half of it) in 
the system of coordinates r, z, the axis « coinciding 
with the axis of symmetry of the nozzle and r being the 
variable radins of the inner surface of the nozzle. The 
direction of the Mach lines is laid off from this curve* 

The direction of the Mach lines is determined from 
the Prandt 1 Busemann method with the aid of a rotating 
semi— ellipse made of a transparent material (celluloid or 
tracing paper). The semiaxes of the ellipse (fig. 30) 

are a = 1 and b = y/fe±i whore k is the specific heat 

V k-1 

ratio (for diatomic gases k = ^p/^v ~ 1.405). For the 

value a = 1 any linear magnitude (100, 200 mm, etc.) may 
be assumed with the condition, however, that the scale 
of the ellipse and velocity diagrams should be the same. 
The scale of the geometric parameters of the nozzle is not 
connected with the scales of the ellipse and velocity dia- 
grams. In general the entire construction, both the vector 
diagram and the axial cross— sect i on of the nozzle may be 
made in absolute units. 

In constructing the wall of the fourth part of the 
nozzle the condition is assumed that the streamlines at 
the initial section of this part (section x) are straight 
and radially directed from the point of intersection of 
the nozzle axis with the direction of the walls of the 
third part. The potential line in this case will be an 
arc of a circle with its center at the point of intersection 
of the streamlines. The radius of this arc in the coordi-, 
nates VpVz will, therefore, be equal to the chosen ratio 
^x^^c (assumed scale) and in the coordinates r, z its 
value is determined by the diameter of the critical section 
and the divergence angle of the nozzle. 

After determination of the magnitude and direction of 
the velocity in the intial section of the fourth part of 
the nozzle, proceed to the graphical construction of the 
walls of this part. 

Tot greater clearness the illustrative computation 
with numerical data is made. The nozzle is designed to 
give a flow with the following parameters^ 
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Pa = 


1.033 kg/cm 


static pregsuro in the flow 


Ta = 


288^ abs. 


temperature of the air in the flow 


a 


0.815 m^/kg 


specific volume of the air in flow 


w = 


900 ml sec 


velocity of the flow 


G = 


1^5 kg/scc 


rate of discharge of air through the 
nozzle 


Then 


f = _a 
^ w 


1.5 X 0.815 ^ 
= -JL. = 1,36 X 10^ m^ 

900 



whence 

= 41.6 mm ^ 42 mm ^ 



The. thermodynamic parameters of the air in the chamber 
(before expansion in the nozzle) is determined 

T o 
T. = — ^- = 691 abs. 



1 

v.=v c'^ = 0^09125 m^/kg 
• a 

.2 -^3 5 

p =:j(^Z^^ +1| =21.5 atm. 

i LV76 0; J 

where 

C = la and k - 1.401 
Pi 

Since w^ is small by comparison ^ith and w^^ 

it is neglected In vrhat follows. 

Knowing the initial conditions the geometric and 
thermodynamic parameters in the critical section Tith the 
aid of formulas (1) to (5) may be determined. Carrying out 
the corresponding computations 

f^ = 0.000454 m^ 
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hence 



0.785 



= 24 mm; 



w„ = 481 m/sec 
c 

= 11.35 kg/ca' 
V = 0,144 ffl^/kg 

C 



To determine the parameters at the start of the 
fourth part, that is, at section x there is assumed 



= 1.05, i.e. 

c 



yr^ = 1.05 w = 505 m/sec 



Then 



Px = Pi 



1 - 



2g 



k-1 



= 10.6 kg/cn® 



k-1 

1 

v_ = V. f'^\ = 0,1505 m /kg 



^ Vp 




The area at section x will be 



f = 



^ 456 mm^ 



whence 



dy = 24.2 mm 



1. Proceeding to the graphical construction of the 
wall of the fourth part of the nozzle prepare a semiellipse 
of the type shown on figure 30 on a special sheet of tracii 
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paper. I.t is ass\imed_t_bat a 100 milllmetera and 



"b = 245 millimeters. 

2. On a separate sheet of paper draw the coordinate 
axes Vj,, (the vector diagram of the velocities, fig 

31). From point 0 draw a straight line at the angle 



center 0 an arc intersecting the v^ axis, is described. 

The segment of the arc included in the sector with angle 
a/2 gives the velocity distribution at the initial section 
of the fourth part. On this arc points 1, 2, and are 
arbitrarily denoted. The lines 0—1, 0—2, and 0—3 are 
the velocity vectors for this section. 

Simultaneously with the construction of the vector 
velocity diagram construct the upper half of the axial 
section of the nozzle (fig. 32) hearing in mind that the 
points on the graph of figure 32 and the end of the corre- 
sponding vectors in figure 31 are denoted by the same figure^ 

S# Determine the Mach lines passing through the points 
1, 2, and 3 (fig. 32). For this purpose the semielllpse 
(on tracing paper) is placed with its center (point 0) on 
point 0 on figure 31 (fastened with a pin). On rotating 
the semiellipse about its center find the position for 
which the curve of the ellipse falls on point 1 (there 
may be two such postions), the direction of the Mach line 
is drawn parallel to the major semiaxis (b ) through point 
1, figure 32. In the same manner determine the direction 
of the Mach lines through points 2 and 3. At these points 
the Mach lines are drawn for two directions according to 
the two possible positions of the major axis of the semi- 
ellipse. Point 4 (fig. 32) is obtained as the point of 
intersection of the Mach lines through points 2 and 3 and 
point 5 as the point of intersection of the Mach lines 
through points 1 and 2. 

4-# ;The velocities at the points 4 and 5 are determined 
in the following manner: through point 2 (fig. 31) draw a 
straight line parallel to the second Mach line passing 
through the point 3 (but not 3 to 4). On the line drawn 
through point 2 (fig. 31) take a segment equal to 6 = 

a Vyds/w r^ where w is the velocity at point 2 (equal 




uni t s , that is , 



a/2 = 6^. iTith radius 



r = 1.05 = 105 millimeters and 
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to the vector 02), a = 




2^2 



the velocity of 



adnnd at this point , v^, the component of the velocity 

along the normal to the axis, r the distance of point 2 
(fig. 32) from the axis, and ds the magnitude of the 
segment 2^? C^ig* 32)^ For point ? the segment 8 is 
determined in the same way. In computing a the value 
of w is determined from the vector diagram in linear 
units and to the same scale to which the diagrams are 
drawn. The magnitudes r and ds are determined on the 
diagram in millimeters; Vg is measured from the graph 
of figure 31 with account taken of the scale. 

The question as to the direction in which the segment 
6 should "be laid off is solved hy the following; rule of 
signs. The segment 8 is laid off to the right, that is, 
the projection of the vector 6 on the v^ axis has a 
positive direction in case the product v^ d^ (v^^-a^) 

is positive; otherwise the segment 6 is laid off to the 
left. 

The magnitude d^ is the axial component of the 

velocity (fig, 32), Tor example, for point 2 (fig, 32) 
dg is equal to the segment 2—4 (in the general case d^ 

is equal to the vector from the point at which the velocity 
is known to the point where it is sought). For the point 
2 the magnitude is positive, the values v^.^ — a^ and 

dg are also positive and therefore the segment is positive^ 
The same is true for point 3. 

At the ends the segments 6 are drawn perpendicular 
with intersection point 4 (fig. 31). Connec^j.ng point 4 
with the origin of coordinates the vector 0^4 which 

determines the velocity at point 4 (fig. 31) is obtained. 

To determine the velocity vector 0—5 (fig. 31) it is 
is necessary to determine the value of the segment 6 at 
point 1. For the center of the jet r = 0 and Vj, = 0. 

Hence to determine point 5 the values of Vj. and r are 
assumed the same as for the neighhoring point 2. 

5. Determine, hy the method already described, the 
Mach lines through the points 4 and 5 (fig, 32), The 
point 5 (fig. 32) is determined as the point of intersec- 
tion of the Mach line through point 5 with the z axis 
and point 7 as the point of intersection of the Mach lines 
through points 4 and 5, 
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6m Determine the velocities at points 5 and 7. 
Through point 4 (fig^ 31) draw a straight line parallel 
to 47 (fig^ 32) and parallel to the Mach line through 
point 4_(fig»_32)^ Draw . through point 5 (fig^ 31) par— 
alels 57 and 56 (fig* 32). On the straight lines lay off 
the segments 6, Prom the ends of the segments 6 draw 
straight lines normal to them, the -ijitersection of which 
gives points 6 and 7, The vectors 06 and 07 are the vector 
velocities at points 6 and 7 (fig, 32). 

7. Determine the Mach lines through points 6 and 7 
(fig. 32). Point 8 is determined as the point of intersec- 
tion of the Mach lines through points 6 and 7. 

8. Determine the velocity at point 8 (fig, 32)# For 

this purpose draw through point 7 (fig, 31) a straight 

line parallel to the Mach line through point 7 (hut not 78); 
through point 6 (fig, 31) a parallel to the Mach line through 
point 6 (hut not F5), On the straight lines lay off the seg- 
ments 6. It should he noted that at point 6 the value v^ = 

0 = r hence the value of 6 for point 6 is taken correspond- 
ing to point 7. Prom the ends of 6 draw straight lines per- 
pendicular to 6. Their intersection will give point 8 (fig, 
31), 

a. Determine the Mach lines through point 8 (fig. 32). 
Point 9 is determined as the point of intersection of the 
Mach line through point 8 with the axis, 

10. Determine the velocity at point 9, For this pur- 
pose through point 8 (fig, 31) draw a straight line parallel 
to the line through point 8 (fig, 32) (not 89), On this 
straight line lay off the segment 6» Prom the ends of the 
segment 5 and perpendicular to the latter draw a straight 
line whose intersection with the Vg axis will give point 9. 

Having determined point 9 the character of the flow in 
the entire triangle of definition 1, 3, and 9 is obtained, 
that is, that part of the flow Trhich is entirely determined 
by the distribution of the velocity in the initial section 
and on which the part of the 'rail lying to the right of the 
initial section has no effect. 

Hext solve the problem of what shape it is necessary 
to give the nozzle wall for a given velocity distribution 
at the initial section in order to obtain, at the outlet 
section, a velocity constant in magnitude and direction^ 
This velocity should be equal to the velocity at the point 



NACA TU No^ 1066 



31 



9 and, therefore , the Mach line from point 9 will be straight^ 

11, Draw the Mach line through point 9 (fi^; 72)^ On 
it take a number of points denoted by 9 with a subscript^ 
The velocities at these points will all agree #ith the 

vector 09 (fig^ 31)^ 

13. Determine the point 10 which is the point of 
intersection of the Mach line through 8 and the second 
point 9, 

13. Determine the velocity at point 10. To do this, 
draw through point 8 (fig* 31) a straight line parallel to 
89 (figt 32) and on it lay off the segment 6. From the 
end of, the latter draw a perpendicular line to intersect 
with the straight line from point 9. The point of inter- 
section will be the end of the velocity vector at point 
10. 

14. Determine the point 11 (fig. 52), which is the 
point of intersection of the Mach lines through points 
7 and 10. 

15. Determine the velocity at point 11. The vector 
0—11 is the velocity vector at point 11. 

16. Determine point 12 (fig. 32), which is the point 
of intersection of the Mach lines through points 4 and 11. 
Determine the velocity at point 12. For this purpose draw 
through poJ.nt 4 (fig. 31) a straight line parallel to tho 
segment 47 (fig, 32) and through point 11 a line parallel 
to tho segment 7 - 11. The vector 0-12 determines the 
velocity at point 12. 

17. Determine* point 13 (fig* 32), the point of inter- 
section of the Mach lines through points 10 and 9, 

18. Determine the velocity at point 13 (fig. 32)^ 
The vector 0—13 is the velocity vector at point 13 
(points 13 and 10 (fig. 31) practically coincided). 

19. Determine point 14 (fig. 32), which is the point 
of intersection of the Mach lines through points 11 and 13. 

20. Determine the velocity at point 14. For this pur- 
pose through point 11 (fig. 31) draw a line parallel to 
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10 - 11 (fig* ?2) a nd t hrough point 13 (fig, 31) draw 
a line parallel to 10 - 13 (fig* 32). 

21, Determine the point 15 (fig. 32) which is the 
point of intersection of the Mach lines through 13 and 9. 

Thus a sufficient number of points in the angular 
space 399. have been determined. 

22, The fi:enerator of the wall is obtained in the 
following manner: through point 3 (fig. 32) draw a straight 
line in the radial direction, that is, pr olon£_^he_wall of 
the third part to intersect with the segment 4—12 and the' 
point of intersection a is obtained. 

The velocity at point a is_obtained by linear inter- 
polation, that is, the segment 4—12 (fig. 31) is diy ided 
in the same ratio in which a divides the segment 4—12 
(fig, 32). The vector oa (fig. Si) thus obtained gives the 
direction of the followj.n£__ segment of the wall generator to 
the intersection with 11 - 12 Fhich is drawn in figure 32 
parallel to the segment oa , and so forth, prolonged to 
intersect with the line 99 and the latter point of inter- 
section gives us the outlet section. For greater clarity 
the point, a and similar points are not shown on the dia- 
gram. 

The results of the computation by the graphical con- 
struction are collected in a table. The velocity w, ob- 
tained at the outlet section, may be less than or greater 
than the velocity which must be obtained since the ratio 
w^/wq was arbitrarily assumed. Assuming, successively, 
the values w^/w^ = 0.2, 1.5, 1.8, and so forth, w/w^ is 
obtained for each ratio ^x/^c* ?rom the obtained values 
of w/w^ for different w^^/w^ draw the curve of w/w^ as 

a function of w^/w^ (fig. 33), From the curve find the 
values of w^^/w^ required for the given velocity. As may 
be seen from the curve, for w/w^ = 900/481 = 1,87 the 
value of ^x^^c ~ It should be noted that the curve 

on figure 33 is true only for the divergence angle a=12^ 
and for the discharge rate G = 1.5 kg/sec assumed in the 
computation. The vector velocity diagram and the nozzle 
section obtained for the given velocity -rith w2^/w^=l,67 
are shOTrn in figure 34, The equivalent Laval nozzle (the 
straight line generating tho Trail) is shown dotted. 

Translation by S, Reiss, 
National Advisory Committee 
for Aeronautics, 
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Figure Apparatus for measuring 

the impulse of a jet with 
dynamometer and disk* 
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Figure 2*— Sketch of apparatus for measuring the reaction 
force of a jet. 




Tigure 3r 
General 
view of 
the appa- 
ratus 
sketched 
in Fig. 2. 
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- -Section fhroaqh Ad 




Figure 4.- 
A battery 
of Pitot 
tubes. 



Figure 10#- 

Sctieme of 
two-dimen- 
Bional 
supersonic 
flow for 

p > P • 

^a^ ^0 



Figiire 11.- 

Scneme of 
two-^dimen- 
sional 
supersonic 
flow for 

Pa<Po- 
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Figure 16.- 
Sketches of 
nozzles for 
studying the 
effect of the 
inlet part. 



Figs. 5,7,8,9 
12 
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Figure 7.- Velocity distribution 

in the transverse 
cross^section of a jet from a 
conical Laval nozzle. 
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nj/ sec 

Figure 5.- Dynamic pressure in 

supersonic flow com- 
puted by the Baleigh formula. 
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Figure 8.- Velocity and 

static 
pressure distribution 
along the axis of a jet 
from a Laval nozzle. 
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Figure 9.- Two-dimensional supersonic flow. 
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Figure 24.- Flow spectriun for Figure 25.- Flow spectrum for 
the conical nozzle. the Frankl nozzle 



Figs. 13.14.15 . 
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Figure 13.- Sketch of conical 
nozzle, 



Figure 14.- Dis- 
tribution of 
dynamic pressure 
for two nozzles 
with different 
shapes of in- c 
let parts 
(x = 20 mm) . 
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Figure 15.- Dynamic pressure 

and velocity 
distributions for two nozzles 
with different shapes of inlet 
parts (along the flow axis in 
the working part) . 



D = Diameter of flow. 
X - Distance from nozzle 
outlet. 
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Figs. 17,18,20 




Figure 17.- Curves of variation of the static pressure and ve- 
locity for the Stanton nozzle. 

kg 



Figure 18.- Static 

■ pressure and 
velocity distributions 
along the flow axis for 
the conoidal nozzle. The 
nozzle was computed for 
w = 1.82c. 
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Figure 20.- Contour of nozzle constructed 
"by the method of Prof. Frankl. 
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Figs. 22,23 
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7i^re 23.- Transverse velocity 
distrilDution for 

w/c = 2.64. 



Figure 22.- Transverse velocity distribution 
for w/c = 1.47, 
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Figure 26.- Dependence of the reaction 

of the jet on the diver- 
gence angle of the nozzle. 



,'Pth 



20 



p 10 



1\ 


















































[ . 


— d 


y 


1 






— c 


) — 






— < 


) 
























■H 


1 























































1 1.2 1.4 1.6 1.8 2.0 2.2 

de/dc 

?igure 27.- Dependence of the reaction of the 
.jet on da/lcf. 
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Figure 28.- Variation of mw and Ap with p^, 
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Jigs. 29,30,33 
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Pigiire 33.- Curve of dependence of 
w/wq on T/^/^^c• 
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Figs. 31,33,34 
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Figure 34.- Vector velocity diagram and axial section of the 
nozzle for w ta 900 m/seo. 
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